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Abstract 

We investigate correlations between fluorescence photons emitted by 
single N-V centers in diamond with respect to the optical excitation 
power. The autocorrelation function shows clear photon antibunching 
at short times, proving the uniqueness of the emitting center. We also 
^— report on a photon bunching effect, which involves a trapping level. An 
■ analysis using rate-equations for the populations of the N-V center levels 

shows the intensity dependence of the rate equation coefficients. 

a,: 
^— > 

Introduction 

Quantum cryptography relies on the fact that single quantum states 
rN \ can not be cloned. In this way coding information on single photons 

would ensure a secure transmission of encryption keys. First attempts 
for quantum cryptography systems [1, 2] were based on attenuated laser 
pulses. Such sources can actually produce isolated single photons, but 
the poissonian distribution of the photon number does not guarantee 
both the uniqueness of the emitted photon and a high bit-rate. Therefore 
a key milestone for efficient and secure quantum cryptography systems 
is the development of single photon sources. 

Several pioneering experiments have already been realized in order 
to obtain single photon sources. Among these attempts one can cite 
twin-photons experiments [3, 4], coulomb blockade of electrons in quan- 
tum confined heteroj unctions [5], or fluorescence emission from single 
molecules [6, 7]. Up to now interesting results were obtained, but con- 
siderable work is still needed to design a reliable system, working at room 
temperature with a good stability and well-controlled emission proper- 
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ties. More recently several systems have been proven to be possible 
candidates for single photon sources. Antibunching was indeed observed 
in CdSe nanospheres [8], thereby proving the purely quantum nature of 
the light emitted by these sources. In the same way photon antibunching 
in Nitrogen- Vacancy (N-V) colored centers in diamond was reported by 
our group [9] and by Kurtsiefer et al [10]. A remarkable property of these 
centers is that they do not photobleach at room temperature: the fluo- 
rescence level remains unchanged after several hours of continuous laser 
irradiation of a single center in the saturation regime. These centers are 
therefore very promising candidates for single photon sources. 

In this paper we present further investigation of the fluorescence light 
emitted by NV colored centers. We analyze the system dynamic through 
the measurement of the autocorrelation function, showing the existence 
of a shelving effect which reduces the counting rate of the fluorescence 
emission, and gives rise to photon bunching [11]. We study the de- 
pendence of this behavior on the pumping power, and we compare the 
experimental results with the predictions of a three-level model using 
rate equations. 

1. EXPERIMENTAL SETUP 

Our experimental set-up is based on a home-made scanning confocal 
microscope. 




Figure 1 (a) Confocal microscopy raster scan (5x5 //m 2 ) of the sample performed 
about 10 fj,m below the diamond surface. The size of a pixel is 100 nm. The integration 
time per pixel is 32 ms. The laser intensity impinging on the sample is 15 mW . (b) 
Line scan along the dotted line. The data is shown together with a gaussian fit, 
which is used to evaluate the signal and background levels (right). Here we obtain 
p = S/(S + B) =0.81 
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A CW frequency doubled Nd:YAG laser (A = 532nm) is focused on 
the sample by a high numerical aperture (1.3) immersion objective. A 
PZT-mounted mirror allows a x-y scan of the sample, and a fine z-scan 
is obtained using another PZT. 

The sample is a 0.1 x 1.5 x 1.5mm 3 [110] crystal of synthetic lb diamond 
from Drukker International. The Nitrogen- Vacancy centers consist in a 
substitutional nitrogen with an adjacent vacancy, and are found with a 
density of about 1 [im~ 3 in lb diamond. The centers can be seen with a 
signal to background ratio of about 5 by scanning the sample as shown 
in figure 1, and a computer-controlled servo- loop allows to focus on one 
center for hours. 

The fluorescence is collected by the same objective, and separated 
from excitation light by a dichroic mirror. High rejection (10 10 ) high- 
pass filters remove any leftover pump light. Spatial filtering is achieved 
by focusing on a 50fim pinhole. The fluorescence is then analyzed by 
an ordinary Hanbury-Brown and Twiss set-up. The time delay between 
the two photons is converted by a time-to-amplitude converter (TAC) 
into a voltage which is digitalized by a computer board. 




Figure 2 Experimental data of g '(t) for different values from 0.3 to 31 mW of 
pumping power. The solid line represents the fit. The dashed (resp. thin) lines are 
calculations of g (2 \r) for &23 and £32 (see figure 3 for definition) not depending (resp. 
depending) on pump power. Details are given in the text of section 2. 
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In a low counting regime, this system allows us to record the autocor- 
relation function: g^{r) = <^gHg±£> . 

Experimental data were corrected from background noise, and nor- 
malized to the coincidence number corresponding to a poissonian source 
of equivalent power [9]. Different correlation functions obtained for sev- 
eral pumping power are shown in figure 2. These functions are centered 
around r = 0, but one can notice slight drifts attributed to thermal drift 
of our acquisition electronics. For r = 0, g( 2 \r) goes clearly down to 
zero, which is the signature of a single emitting dipole. On the other 
hand, g^ 2 \r) goes beyond 1 at longer times, and then decays to 1. This 
behavior is an evidence of the presence of a trapping level. To correctly 
describe the system we therefore considered a 3-level system. 

2. THEORETICAL BACKGROUND AND 
DISCUSSION 

Owing to the fast damping of coherences, we use rate equations. Let 
us consider the 3-level scheme described in figure 3. The evolution of 
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Figure 3 Three level system as used for modelisation 
the populations are therefore given by: 
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with the initial condition a± = 1, a% = as = at t = 0, which means 
that a photon has just been emitted and the system is therefore prepared 
in its ground state. The decay rate from level 3 to 1 is neglected [12]. 
By analytically solving equations 1.1, one can derive 

(2)/ \ °2(t) , 1+Se k tm +k lm 1 - g e _Hm-ki m 

9 ( >(t) = 7—7 = 1 — exp 2 _ exp 2 (1.2) 

£72 (OO J 2 2 

with the stationary population: 

( \ ^12^32 n q\ 

K12K23 + K12K32 + hlh2 
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and with 

hm = ku + k 2 i + k 2 3 + k 32 (1.4) 
kf m = (( fc i2 + k 2 i - k 23 - k 32 ) 2 + 4fe 2 ife 23 ) (1.5) 

+ k 32 (ku + k 21 - k 23 - k 32 ) 

9e = u — L ( 1 - 6 ) 

k\ m k 32 

We thus have four equations which enables us to express the four rates fey 
with respect to the experimental variables g e , kt m , ki m ,a 2oo . By fitting 
the experimental values of with expression 1.2 we obtain the values 
of g e , k t m and k\ m for each value of the intensity. The last value needed 
in order to solve the system concerns the parameter a 2oo which is directly 
linked to the count rate 



N = 7] X k 21 x G 2oa (1.7) 

We assume that k 2 \ doesn't depend on the pump power, and we set rj 
to the value for which this condition is satisfied. We found rj = 3 x 
10 -3 which is in perfect agreement with our estimated value [9]. The 
parameters fey are plotted as functions of pumping power in figure 4. 



0.4- 






• Pump */ 




■ k.21 


0.3- 


M k23 // 




A k32 


C 0.2- 


• / 


0.1- 




o.o- 






10 20 30 




mW 



Figure, 4 Evolution of kij with respect to the pump power 



Note that k± 2 is, as expected, a linear function of the pumping power, 
and that fe^ 1 has a constant value of 11.6ns, which corresponds to results 
reported in literature [13]. 

What is noteworthy here is that k 23 and k 32 linearly depend on pump 
power, with k 23 greater than k 32 , which means that at high pumping 
power the system tends to be shelved into the third level. This intensity- 
dependent effect has not been reported yet. In Figure 2, three plots were 
superimposed to the experimental results. The solid line represents the 
best fit. The dashed line represents the result of equation 1.2 with 
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values of ^23 and independent on intensity taken from reference [10]. 
One can see in this last case that the agreement between calculated 
and experimental value is not fully satisfactory. Finally, the thin line 
is obtained by using the values of ^23 and ^32 given by the linear fit of 
figure 4. Agreement with the experimental results is good. 

This intensity dependence of £23 and ^32 can also be observed on the 
total photon counts. In fact the number of photons emitted per second 
should decrease as the trapping in the metastable state increases. Figure 
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Figure 5 Saturation behaviour of the N-V center. The solid line gives the fit by 
using our model, and the dashed line the results for normal saturation behaviour. 
Values are plotted in respect to pumping power. The intensity may be deduced from 
the focused beam size, which in our case is not accurately known due to spherical 
aberrations. 

5 shows the count rate of the single N-V center as a function of the pump 
power. We clearly see a decrease in the photon counts. We can fit with 
a good agreement our experimental data using equations 1.7 and 1.3. 

We have repeated this experiment with the 532nm Nd:YAG line and 
the 514nm Argon line and found similar linear dependency of the ^23 
and &32 rates with the pump power. 

Conclusion 

We have measured the autocorrelation function of the fluorescence 
light emitted from a single NV colored center in diamond using a confocal 
microscope. Our results are in very good agreement with what can be 
expected from rate equations in a 3-level scheme, provided that a linear 
dependence of some rate coefficients with pumping power is assumed. 
Such an effect has not been reported yet, and must be taken into account 
when designing a pulsed single-photon source. 
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